Screening in two-dimensional foams 
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Using the Surface Evolver software, we perform 
numerical simulations of point-like deformations 
in a two-dimensional foam. We study perturba- 
tions which are infinitesimal or finite, isotropic 
or anisotropic, and we either conserve or do 
not conserve the number of bubbles. We mea- 
sure the displacement fields around the pertur- 
bation. Changes in pressure decrease exponen- 
tially with the distance to perturbation, indi- 
cating a screening over a few bubble diameters. 



1 Introduction 

A foam is a discrete material, made of gas bub- 
bles separated by a continuous liquid phase. 
Foams act as elastic solids for small deforma- 
tions, but when large strains are applied they 
behave as viscous liquids; applied stresses are 
relaxed by discrete rearrangement events that 
occur in the foam. Other changes are due to 
ageing, where some bubbles gain gas at the ex- 
pense of others. In most cases of foam evolu- 
tion, continuous changes, such as infinitesimal 
changes of bubble shapes and sizes, alternate 
with discontinous processes [I]. 
For instance, the length of an edge might de- 
crease (or, inversely, it might stretch). If it van- 



ishes, a neighbour-swapping event occurs, and 
a new edge is created: this is the topological 
Tl process [2J (its inverse is a Tl too). Alter- 
natively, an edge breakage leads to bubble co- 
alescence, also called fusion (its inverse is a di- 
vision). All these perturbations are anisotropic. 
On the other hand, when a bubble's area de- 
creases, it is an isotropic pertubation (its in- 
verse is a bubble inflation). If it vanishes, a 
reduction in the number of bubbles occurs: this 
is the topological process T2 [2J, reminiscent of 
biological cell apoptosis (its inverse is a nucle- 
ation). 

Each pertubation affects the neighbouring bub- 
bles, over a certain range. This range has been 
the subject of various studies, with diverse mo- 
tivations, all in two dimensions. The perturba- 
tion induced by a Tl has been measured in ex- 
periments [1] and in simulations [5]. The effect 
of changing the volume of a single bubble is also 
studied, both in experiment [6j [7] and in sim- 
ulations [HI [9]. A laser has been used to break 
the wall between bubbles [ID] . This laser abla- 
tion is also used in biological tissues [IT] , and 
cell division and apoptosis are currently studied 
too [12] , motivated by a desire to understand 
the mechanics of biological tissue. 

An infinitesimal change is expected to have 
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an infinitesimal effect on neighbouring bubbles. 
Conversely, finite discontinuities can have a fi- 
nite effect on other bubbles. If these events 
are scarce enough, and spatially screened by 
the foam disorder, their effects on the foam are 
independent, and the foam can at large scale 
be treated as a continuous material [13] . The 
foam properties of insulation against sound and 
explosion probably derive from this screening. 
A question under debate is to determine under 
which conditions these processes correlate, trig- 
gering for instance a cascade of Tls [H HI [7J [9] . 
This can lead to large-scale fluctuations such 
as avalanches [HI [15], [161 E] • The role of dis- 
order on perturbation screening is also of po- 
tential interest in various other discrete systems 

ps USES. 

To study the screening in foams, we focus on a 
2D simulated foam with an isolated perturba- 
tion. Physically, this is relevant for a system 
which is initially formed by well separated de- 
fects, that may eventually interact. We try to 
distinguish a screening area around the pertur- 
bation, and the rest of the foam where no dis- 
turbance is felt. 

We simulate small continuous geometrical 
changes, then finite discontinuous topological 
processes. We either conserve or change the 
number of bubbles. We try different foam dis- 
orders and boundary conditions. We perform 
isotropic (scalar) pertubations, and anisotropic 
(tensorial) ones. Two types of measurements 
are extracted from the simulations: the dis- 
placement of each bubble (a vector), and the 
change in pressure [21] in each bubble, with re- 
spect to the inital (unperturbed) foam. 

2 Method 

We create a disordered foam structure using 
a Voronoi construction based upon a random 
Poisson process for generating the cell centres 
[22] ■ The cell areas A are chosen based upon a 
randomly generated Weibull distribution 
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1.115(A), so that the peak of the distribution is 
close to A = (A). For the monodisperse case, 
each cell is set to have the same target area, 
rather than taking the limit (3 — > oo. We also 
constructed ordered monodisperse foam consist- 
ing of 8 and 20 concentric shells of hexagons [21] . 
The disorder of each foam is defined by the (di- 
mensionless) variance of the area distribution: 



var(A) 



(A-(A)f 
(A) 2 



(2) 



where the parameter j3 > 1 determines the area 
dispersity and the parameter A is chosen as A = 



where () denotes an average over the whole 
foam. 

Each foam is first equilibrated in the Surface 
Evolver software [25], using a mode in which 
each film is represented exactly circular arc 
and a value of surface tension 7 equal to one. Tl 
topological changes are automatically triggered 
when films shrink to a very small length during 
the equilibration process. Some of the initial 
foam structures used are shown in figure [U 
We then trigger changes to the structure as fol- 
lows, illustrated in figure [21 

• An edge stretch is performed on the longest 
edge of the bubble nearest the centre of the 
cluster, by extending or reducing its length 
by a fraction 5L and fixing its endpoints. 

• An inflation of the bubble nearest the cen- 
tre of the cluster by a fraction 5 A. 

• A Tl is performed on the shortest edge of 
the bubble nearest the centre of the cluster. 

• A T2 event is triggered by removing the 
bubble nearest the centre of the cluster. 

• A coalescence event is triggered by remov- 
ing the longest edge of the bubble nearest 
the centre of the cluster. 

• A division event is triggered by inserting a 
edge across the centre of the bubble nearest 
the centre of the cluster, giving two bubbles 
each with half the area. 

In each case we re-equilibrate the structure; we 
allow Tls to occur, but in fact they occur in- 
frequently, except after the T2. We record the 
position of the bubble centres (average of the 
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vertex positions) and the bubble pressures be- 
fore and after the distortion, denoted Xi, y^ pi 
with i — 0, 1. 



3 Bubble displacement 

For free clusters, the displacements of the bub- 
ble centres after each event are shown in fig- 
ure [31 Since we apply no constraint during 
the equilibration to prevent rigid body displace- 
ments and rotations, we first calculate the av- 
erage over all bubble centres of the x and y dis- 
placement and the average rotation about the 
centre of the cluster. These averages are then 
subtracted from the displacement of each bub- 
ble. The results suggest that Tls, T2s, coales- 
cence, division and edge- stretching events all in- 
duce approximately quadrupolar displacement 
fields, while the inflation event is almost purely 
radial. It is surprising to see that the effect of 
the T2 is not isotropic, although not surpris- 
ing to see that it has the largest effect (greater 
displacements) . 

This information is summarised in figure HI 
which shows the variation in position of each 
bubble centre as a function of radial and an- 
gular position. The radial displacement de- 
creases with radial distance, and, apart from 
the inflation event, varies sinusoidally with an- 
gle like a quadrupole, Ar ~ sin(26 l ). The an- 
gular displacement field is slightly asymmetric 
but also approximately sinusoidal and dipolar, 
tqA6 ~ sin(#). The exception is the inflation 
event, for which the motion is radially outwards. 

This dipolar angular displacement is not 
screened and is sensitive to boundary condi- 
tions. In foams with periodic boundary condi- 
tions the angular displacement deviates from a 
sinusoid (figure E]) . In foams with fixed bound- 
aries (figure [H bottom), which force the radial 
displacement to vanish (equivalently, each per- 
tubation interacts with its reflection at bound- 
aries, or virtual image), the angular displace- 
ment becomes quadrupolar. 



4 Pressure 

The perturbation induces a change in each bub- 
ble pressure, Ap = p 1 — p , with respect to the 
initial foam. We normalize it by 7/y (A) to fa- 
cilitate comparison between different foams. It 
decreases with radial distance r = \Jxq + ?/q, 
as shown in figure EJ 

A small change in the area of a bubble at the 
centre of the cluster of 150 bubbles, in this case 
by a factor of 20%, has only a small effect on 
the pressure differences, which decrease roughly 
exponentially with radial distance. A similar 
result was found in a cluster of 1400 bubbles. 
All other distortions induce a stronger variation 
in pressure difference. The largest pressure dif- 
ferences are associated with the T2 and coales- 
cence events. In the case of a change in edge 
length, there is little difference between exten- 
sion and compression; similarly, inflating and 
shrinking a bubble have the same effect. 
The pressure change fluctuates from bubble to 
bubble, and is either positive or negative. Its 
standard deviation visibly decreases (figure E]). 
In an ordered cluster, it is possible to define ra- 
dially concentric shells of bubbles and examine 
the pressure change in each cell, but to quan- 
tify the pressure differences in disordered foams 
the shells vary greatly in composition [3] and 
it is necessary to bin the data in some manner. 
We cho ose a b in width of one bubble diameter 
D = 2^J (A)/tt, and increase the position of the 
centre of the bin in intervals of -D/10. The aver- 
age pressure difference Ap is calculated in each 
bin: it is always close to zero beyond the central 
bubble. The normalised standard deviation in 
each bin, std(Ap), directly compares the effect 
of the different pertubations we apply (figure [7]). 
Remarkably, all pertubations decrease with a 
similar exponential decay, sometimes over more 
than two decades. Its characteristic length, or 
screening length £ s , is the inverse of the slope: 

1 _ dlog(std{Ap)) 

i s ~ dr ■ {6) 

Although we did not find quantitative measure- 
ments in the literature, the available qualitative 
data [H [TTJ Q2] seem to observe a screening over 
one or a few bubble diameters. This is what 
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we observe for the different distortions (figure 
[8]): despite large changes in several parameters, 
there is little variation in £ s , with no clear de- 
pendence upon cluster size (data not shown) 
or disorder. It is in general larger for confined 
foams and those with periodic boundary condi- 
tions. 

The choice of var(A) as the disorder parameter 
requires some discussion: since the screening is 
a local effect, we should ideally use a local mea- 
sure of disorder. It is, however, not clear how 
to define such a measure. 

For all distortions except the bubble inflation, 
the exponential decrease of std(Ap) crosses over 
to a plateau near the outer part of the foam. 
The cross-over length r c marks the limit where 
the effect of the pertubation can be detected. 
We measure it as the intersection of linear fits 
to the decrease and the plateau (figure [7b). Its 
value is very robust (figure Mj: it is roughly 60% 
of the cluster radius independent of pertuba- 
tion type and boundary conditions; as well as of 
polydispersity and cluster size. In the present 
simulations, it corresponds to roughly 5 to 15 
bubble diameters. Thus, we can change the 
volumes of bubbles around the remaining 40% 
outer part of the foam: it does not change the 
screening length (figure [7b). 



4.1 Other scalar measures of 
screening 

The perturbation induces other changes in 
the structure of the foam, for example each 
edge may shrink or lengthen to accommodate 
changes at the centre of the cluster. We jus- 
tify a posteriori our choice of pressure difference 
as a measure of screening by comparing it with 

(i) the change in each bubble's perimeter, Ae, 
normalized by the square-root of its area, and 

(ii) the change in the length of each edge A£ 
[26] , normalized by <J (A). Figure [10] shows that 
neither of these scalar quantities indicates any 
screening in both ordered and disordered foams, 
while the pressure change does. In particular, 
A£ ~ r~ 2 over more than a decade. 



5 Summary and perspec- 
tives 

In simulated 2D foams, we have performed con- 
tinuous and discontinuous perturbations, either 
isotropic or anisotropic, which conserve or not 
the number of bubbles. We have varied poly- 
dispersity over a decade and cluster size over 
almost a decade, with free, periodic or fixed 
boundary conditions. 

The bubble displacements are dipolar in all 
cases, except for bubble inflation. They extend 
to the foam boundary and are thus sensitive to 
boundary conditions. 

The pressure change in each bubble fluctuates 
from one bubble to the other, even for bubbles 
the same distance from the perturbation. The 
standard deviation within each of these "shells" 
decreases exponentially with distance from the 
perturbation over up to two decades, with a 
characteristic (screening) length of the order of 
one average bubble diameter. 
At a distance close to 60 % of the foam ra- 
dius, the standard deviation of pressure differ- 
ence reaches a plateau. Beyond this distance 
the pertubation does not significantly affect the 
pressures. 

Perspectives include the measurement of ten- 
sorial quantities such as the deformation; the 
comparison with displacement and deforma- 
tion fields predicted by standard elasticity the- 
ory; direct comparison with experiments; vecto- 
rial anisotropic pertubations, such as obtained 
by moving a point-like defect within the foam 
(Stokes experiment [13]); and transition from 
low disorder to the perfectly ordered case (fig- 
ure dH). 
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Fig. 1 Representative images of the foams used. Top row: Monodisperse foams, both ordered 



(217 and 1261 bubbles) and disordered (150 bubbles). Middle row: Polydisperse foams, 
in each case with almost the same value of polydispersity, var(A) pa 0.12, and 250, 450 
and 1400 bubbles respectively. Bottom row: Polydisperse foams without boundaries 
(periodic, 400 bubbles) and with boundaries (383 bubbles) respectively. 



Fig. 2 The six types of distortion applied to the foams, illustrated here for the disordered 



monodisperse foam of 150 bubbles. Edge stretch (by 20%), inflation (by 20%), Tl, 
T2, coalescence, division. 



Fig. 3 The displacement field after (a) 10% inflation, (b) 20% edge stretch, (c) Tl, (d) T2, 
(e) coalescence, and (f) division, in a cluster of 1400 bubbles with var(A) pa 0.12. The 
average displacement and rotation has been subtracted in each case, and the vector 
length multiplied by a factor of 200 (in a), 100 (in b) or 30 (in c,d,e,f). 

Fig. 4 Bubble displacement, for the same data shown in figure [3) (a) The magnitude of the 



radial displacement as a function of radial position, (b) The radial displacement as a 
function of angular position, (c) The angular displacement as a function of angular 
position. 



Fig. 5 Effect of boundary conditions. Bubble displacement field after a Tl in a periodic foam 



of 400 bubbles (a) and in a confined foam of 383 bubbles (b). The average displacement 
and rotation have been subtracted, and the vector length multiplied by a factor of 30. 
Angular bubble displacements in the same periodic (c) and confined (d) foams. 



Fig. 6 The effect on the bubble pressures of each of the distortions in figure [2] applied to a 



disordered monodisperse cluster of 150 bubbles. On linear axes the effect of (a) an 
infinitesimal perturbation is much less than after (b) the discontinuous processes (note 
the different scales). 



Fig. 7 (a) Semi-log plot of the standard deviation of pressure change (see text) against radial 



position measured in units of cluster radius. Same data as in figure [61 The drop-off in 
the data at large distance is caused by inaccuracy in defining the radius of the periphery 
of the cluster, (b) Same plot for two highly polydisperse clusters of 250 bubbles. The 
dashed lines are fits to a piecewise linear function. The difference between the clusters 
is that the bubble volumes beyond eight bubble diameters from the centre have been 
randomly varied, as shown in the inset images of the clusters. This doesn't make a 
significant difference to the screening length. 



Fig. 8 Screening length t s of std(Ap) (eq. [3]), in units of bubble diameter D = 2^J{A)/n , 
after (i) inflation (ii) an edge stretch, (iii) a Tl, (iv) a T2, (v) coalescence, and (vi) 
division. 

Fig. 9 Cross-over radius r c , in units of foam radius, after (i) inflation, (ii) an edge stretch, 
(iii) a Tl, (iv) a T2, (v) coalescence, and (vi) division. It is rarely possible to measure 
it for bubble inflation. 

Fig. 10 Log-log plot of the standard deviation of pressure difference Ap, perimeter change 



Ae and edge-length differences Al against radial position measured in units of cluster 
radius, comparing data for a Tl in a disordered cluster of 1400 bubbles with a Tl in an 
ordered cluster of 1261 bubbles. In both cases, A£ decreases quadratically (the solid 
line has slope -2) over the whole foam. Ae does the same in the disordered cluster, 
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but the ordered case shows a more rapid drop, then little change over the outer part of 
the cluster. The only scalar measure of disorder that shows screening in both ordered 
and disordered foams is the pressure change. 

Fig. 11 The effect on the bubble pressures of a 20% inflation event in a monodisperse ordered 
cluster of 1261 bubbles. In contrast to other perturbations to a monodisperse foam, 
and inflation events in disordered foams, there is no pressure screening: the pressure 
difference increases slightly towards the outside of the cluster. Note the scale, orders 
of magnitude smaller than figure 
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Figure 1: Representative images of the foams used. Top row: Monodisperse foams, both ordered 
(217 and 1261 bubbles) and disordered (150 bubbles). Middle row: Polydisperse foams, in each 
case with almost the same value of polydispersity, var(A) rs 0.12, and 250, 450 and 1400 bubbles 
respectively. Bottom row: Polydisperse foams without boundaries (periodic, 400 bubbles) and 
with boundaries (383 bubbles) respectively. 
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Figure 2: The six types of distortion applied to the foams, illustrated here for the disordered 
monodisperse foam of 150 bubbles. Edge stretch (by 20%), inflation (by 20%), Tl, T2, coales- 
cence, division. 
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Figure 3: The displacement field after (a) 10% inflation, (b) 20% edge stretch, (c) Tl, (d) T2, 
(e) coalescence, and (f) division, in a cluster of 1400 bubbles with var(A) « 0.12. The average 
displacement and rotation has been subtracted in each case, and the vector length multiplied by 
a factor of 200 (in a), 100 (in b) or 30 (in c,d,e,f). 
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Figure 4: Bubble displacement, for the same data shown in figure [31 (a) The magnitude of the 
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after (i) inflation (ii) an edge stretch, (iii) a Tl, (iv) a T2, (v) coalescence, and (vi) division. 
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Figure 9: Cross-over radius r c , in units of foam radius, after (i) inflation, (ii) an edge stretch, (iii) 
a Tl, (iv) a T2, (v) coalescence, and (vi) division. It is rarely possible to measure it for bubble 
inflation. 
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Figure 10: Log-log plot of the standard deviation of pressure difference Ap, perimeter change 
Ae and edge-length differences A£ against radial position measured in units of cluster radius, 
comparing data for a Tl in a disordered cluster of 1400 bubbles with a Tl in an ordered cluster 
of 1261 bubbles. In both cases, A£ decreases quadratically (the solid line has slope -2) over the 
whole foam. Ae does the same in the disordered cluster, but the ordered case shows a more rapid 
drop, then little change over the outer part of the cluster. The only scalar measure of disorder 
that shows screening in both ordered and disordered foams is the pressure change. 
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Figure 11: The effect on the bubble pressures of a 20% inflation event in a monodisperse ordered 
cluster of 1261 bubbles. In contrast to other perturbations to a monodisperse foam, and inflation 
events in disordered foams, there is no pressure screening: the pressure difference increases slightly 
towards the outside of the cluster. Note the scale, orders of magnitude smaller than figure [7J 
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